Theory of non-equilibrium transport in the SU(N) Kondo regime 
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Using a Fermi liquid approach, we provide a comprehensive treatment of the current and current 
noise through a quantum dot whose low-energy behaviour corresponds to an SU(A'') Kondo model, 
focusing on the case A'' = 4 relevant to carbon nanotube dots. We show that for general N , one needs 
to consider the effects of higher-order Fermi liquid corrections even to describe low-voltage current 
and noise. We also show that the noise exhibits complex behaviour due to the interplay between 
coherent shot noise, and noise arising from interaction-induced scattering events. We also treat 
various imperfections relevant to experiments, such as the effects of asymmetric dot-lead couplings. 
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I. INTRODUCTION 

The Kondo effect has long served as a paradigm in the 
field of strongly correlated electron physics. It is perhaps 
the simplest example of a system where many-body in- 
teractions can give rise to highly non-trivial behavior: its 
essence involves nothing more than a localized magnetic 
impurity which is exchange coupled to conduction elec- 
trons in a metal. Despite having been studied for over 40 
years, interest in Kondo physics shows no sign of abating. 
A large part of this continued interest has been fueled by 
recent advances allowing the controllable realization of 
unusual Kondo effects in nanostructures. These include 
multi-channel Kondo effects where there are many 
conserved flavours of conduction electrons: such systems 
can give rise to non- Fermi liquid physics, and have re- 
cently been realized using semiconductor quantum dots 
Another class of exotic Kondo effects are so called 
SU(A^) Kondo effects, where N > 2. Such systems in- 
volve only a single channel of conduction electrons, but 
the effective spin of the impurity and conduction elec- 
trons is greater than 1/2. While such systems are still 
described at low-energies by a Fermi liquid fixed point, 
the properties of this Fermi liquid are modified in sev- 
eral interesting ways compared to the spin- 1/2 case Q. 
The case TV = 4 has receivedparticular attention due 
to its realizability in double 0, d, Q and triple quan- 
tum dots [3 as well as carbon nanotube quantum dots 

Research on Kondo physics has also been spurred by 
the possibility of studying experimentally its behaviour 
when driven out-of-equilibrium, where non-equilibrium 
is either achieved by the application of a drain-source 
voltage across a quantum dot [l^, [3l , or by externally 
radiating a quantum dot p^ . The non-equilibrium in- 
duced by a voltage has been the subject of a number of 
recent theoretical works [H [13, [H, M, 

In this paper, we will focus on a topic which com- 
bines two of the above avenues of Kondo research: we 
will study non-equilibrium charge transport through a 



voltage-biased quantum dot exhibiting an SU(A^) Kondo 
effect, focusing on the low-temperature regime where the 
physics is described by an effective Fermi liquid theory. 
We present calculations for both the non-linear conduc- 
tance as well as for the current noise. As has been 
stressed in a number of recent papers [H, [H, H^, the 
fluctuations of current through a Kondo quantum dot are 
extremely sensitive to the two-particle interactions asso- 
ciated with the underlying Fermi liquid theory. This was 
first discussed in the case of the standard SU(2) Kondo 
effect by Sela et al. , and was even measured for this 
system in a recent experiment by Zarchin et al. [2^. As 
discussed in Refs. [2J, [2^, the situation becomes even 
more interesting for A'^ > 2, as now one must deal with 
the interplay between coherent partition noise (due to the 
zero-energy transmission coefficient through the dot not 
being one) and the interaction- induced scattering events. 
Of particular interest is the case = 4, which can be re- 
alized in carbon nanotube quantum dots. Very recently, 
current noise in such a system has been measured ex- 
perimentally by Delattre et al. [23|, though not in the 
low-temperature Fermi liquid regime we describe here. 

The results presented here both clarify and extend 
those presented in Refs. [23, [2^ , as well as provide de- 
tails underlying the calculational approach. Particular 
attention is given to the role of higher-order Fermi liquid 
corrections, something that was not correctly treated in 
previous works (see erratum, Ref. [11]). We show clearly 
how in the A'' = 4 case, such corrections lead to an ef- 
fective shift of the Kondo resonance with applied bias 
voltage. As a result, the non-linear conductance does 
not increase with voltage, as would be expected from a 
simple picture of the Kondo resonance as a resonant level 
sitting above the Fermi energy. These Fermi-liquid en- 
ergy shifts are absent in the usual N ^ 2 Kondo effect. 
We also describe the experimentally-relevant case where 
there is an asymmetry in the coupling between the quan- 
tum dot and the source and drain electrodes. Such an 
asymmetry has not been investigated thoroughly in pre- 
vious works. 
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The remainder of this paper is structured as follows. In 
Secini we outline the basic description of our model and 
the Fermi liquid approach. Sec. IIIII and IIVI are devoted 
to providing a detailed discussion of our results for both 
the conductance and the shot noise, as well as details 
on their derivation. In Sec. |Vl we summarize our main 
results for the conductance and shot noise of a SU(N) 
Kondo quantum dot, and conclude. 



II. MODEL DESCRIPTION 

A. Kondo Hamiltonian 

We give here a compact synopsis of the quantum dot 
model wc study, and how it gives rise to Kondo physics. 
The dot connected to the leads is described by the fol- 
lowing Anderson Hamiltonian [2^ 

H = Hb + Hi, + Ht = ed'^n„ + U ^ rio-ricr' 

k.a 



H = T,k.a £k{bl„bka + ala^^ka) + Hk whcrc 



k.a 



CL/R,fccr is the annihilation operator for an electron of spin 
a = 1 . . . N and energy Sk = hvpk (measured from the 
Fermi energy ep) confined on the left/right lead, is 
the electron operator of the dot and Ua = d\dcy the corre- 
sponding density. U denotes the charging energy, Sd the 
single particle energy on the dot and ^l/r the tunneling 
matrix elements from the dot to the left/right lead. The 
general case of asymmetric leads contacts is parametrized 
by <L = icos6', <r = tsiwO with 6 = [0,7r/2]. 9 = 7r/4 
recovers the symmetric case. The rotation in the basis of 
leads electrons 



bka 



cos sm 
sin 6 — cos ( 



CL.fco- 



(2) 



decouples the aka operators from the dot variables. The 
Kondo screening then involves only the b^a variables. In 
the symmetric case, = 7r/4, bka and Uka represent 
respectively even and odd wavefunctions with the dot 
placed at a; = 0. 

We consider in this work the Kondo limit where the 
charging energy U is by far the largest energy scale. 
Below this energy, the charge degree of freedom on 
the dot is quenched to an integer value and does not 
fluctuate. For = U{1 — m ~ m/N), the number 
of electrons is ^^.n^ = m. The virtual occupation 
of other charge states by exchange tunneling with the 
leads is accounted for by the standard Schrieffcr- Wolff 
transformation [sot (or second order perturbation the- 
ory). It transforms Eq. ^ to the Kondo Hamiltonian 



Hk — Jk 

ka.k' ,a 



ka 5 



(3) 



and Jk = %j — tM r. This Hamiltonian acts in the sub- 

space constrained by = m. In this paper, we con- 

centrate on the choice = U (1 — m—m/N) for which po- 
tential scattering terms vanish after the Schrieffer- Wolff 
transformation. Including potential scattering in the for- 
malism is possible, for example along the line of Ref. [sij . 
It however remains outside the scope of this work where 
we focus on the asymmetric dot-lead couplings. 

The lead electrons transform under the fundamental 
representation of SU(N). With exactly m electrons, the 
localized spin on the dot transforms as a representation of 
SU(N) corresponding to a single column Young tableau 
of m boxes. A basis of generators for this SU(N) repre- 
sentation is formed by the N'^ — 1 traceless components 
S"^"' = d\d„, - {mlN)5„.„: with {a, a') ^ {N,N). This 
basis can be used [32| to rewrite Eq. ([3]) as an antiferro- 
magnetic coupling 



Hk — JkS ■ T, 



(4) 



between the impurity (dot) spin S = {S^,A = 
1, . . . , N'^ — 1) and the spin operator of the lead elec- 
trons taken at X = 0, f = iT,k,k' ^a,a' ^lata.a'bk'a' , A = 

1, . . . , — !)■ The N X N matrices are generators of 
the fundamental representation of SU(N), while are 
matrices acting on states with m 



Ni. 



ml. {N~m)\ 



m! {N~m)\ 

electrons. 

Starting from high energies, Jk grows under renor- 
malization. It presages the complete screening of the dot 
spin by the formation of a many-body SU(N) singlet in 
the ground state. A large body of studies has shown 
that the strong coupling fixed point that dominates at 
low energy is a Fermi liquid one. Exact results from the 
Bethe-Ansatz [s^ find low energy exponents that charac- 
terize a Fermi liquid. Writing the Kondo Hamiltonian ^ 
in terms of current, Affleck [s^l has shown by complet- 
ing the square that the impurity spin can be absorbed 
by lead electrons. The resulting (conformal field) theory 
is that of free fermions and it is believed to describe the 
strong coupling fixed point. It shows a simple translation 
of energies in the spectrum corresponding to an electron 
phase shift imposed by the Friedcl sum rule 



N 



(5) 



The identification of the leading irrelevant operator at 
this fixed point yields Fermi liquid behavior |3l| . Alter- 
natively and following Ref. [ij the ground state of Eq. (j4]) 
has been shown [s^ to be a singlet state. Turning the 
coupling to the leads does not destabilize this singlet 
leading again to Fermi liquid exponents. Finally, Numeri- 
cal Renormalization Group (NRG) calculations have con- 
firmed this picture for SU(2) [H and SU(4) 
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B. Fermi liquid theory 



reads 



We now discuss in detail the Fermi liquid theory for 
the Kondo effect, first introduced by Nozieres [s^. It de- 
scribes the low energy regime - the vicinity of the strong 
coupling fixed point - and allows one to make quantita- 
tive predictions even in an out-of-equilibrium situation. 
In Ref. 3 , the Fermi liquid theory of Nozieres has been 
extended with the introduction of the next-to-leading or- 
der corrections to the strong coupling fixed point. These 
corrections are necessary in the SU(N) case for observ- 
ables like the current and the noise since their energy 
(ksT, eV or fJ-sB) dependence is mostly quadratic. 

The Kondo many-body singlet (also called the 'Kondo 
cloud') having been formed, we wish to describe how lead 
electrons scatter off it. At low energies, two channels 
open: an elastic and an inelastic one. Both take place at 
the dot position a; = 0. Elastic scattering is described by 
an energy-dependent phase shift. At the Fermi level ep, 
it is equal to 6o, see Eq. ([5]). We expand the phase shift 
around the Fermi energy. 



Tk 



^ K 



(6) 



where the energy e is measured from ep. ai and a2 are 
dimensionless coefficients of order one. 

It is instructive to think of the elastic scattering off 
the Kondo singlet in terms of an effective non-interacting 
resonant level model (RLM), where this effective reso- 
nance represents the many-body Kondo resonance. This 
is the picture of the Kondo effect provided by slave-boson 
mean- field theory [37| , and is an exact description of the 
SU(7V) Kondo effect in the large N limit [38j]. Note that 
for finite N, one must also deal with true two-particle 
scattering off the singlet, something that will never be 
captured by the RLM; we thus only use it to obtain in- 
sight into the elastic scattering properties. In the RLM 
picture, the first two terms in the phase shift in Eq. ^ 
are attributed to a Lorentzian scattering resonance cen- 
tered at Ek = Tk cot Sq with a width cx Tk [1^ . In the 
SU(4) case, one thus finds that the Kondo resonance is 
centered at a distance Sk ~ Tk above the Fermi energy, 
giving a heuristic explanation for the fact that the low- 
energy transmission coefficient through the dot is only 
T = 1/2. The fact that the Kondo resonance sits above 
the Fermi energy is indeed seen in exact NRG calcula- 
tions of the impurity spectral density d, [l3| • 

The low energy expansion of the RLM phase shift 
(5(e) = arctan^j-^^^ also gives the form Eq. ^ with 

cti/cti = cot(5o. Note that there is no apriori reason that 
this relation must hold for the expansion of the true phase 
shift, as the correspondence to a non- interacting resonant 
level is not exact. Despite this caveat, one finds that the 
prediction from the RLM picture is quite good even at a 
quantitative level. The exact relation between ai and a2 
is extracted Q from the Bethe ansatz solution [3^ and 



a2 _ N-2 r(l/7V)tan(7r/iV) 
of ~ N^l 



N J 



cot So , 



(7) 



where i5o is given Eq. In the SU(2) case, or more 
generally for a half-filled dot with m = N/2, a2 = 0, cor- 
responding to a Kondo resonance centered at the Fermi 
level. This is expected for a model where particle-hole 
symmetry is not broken. In the SU(4) case, Eq. ([7]) 
gives dil a2/al ~ 1.11284 instead of 1 in the RLM. As 
expected, the agreement becomes even better at larger TV, 
and RLM result is indeed the N = oo limit of Eq. ([7|) . 

The phase shift in Eq. ^ completely characterizes the 
low-energy elastic scattering off the Kondo singlet. For 
further calculations, it is useful to describe it using a 
Hamiltonian formulation. The free Hamiltonian describ- 
ing purely elastic scattering is given by 



Ho 



(8) 



Decoupled from the outset, the aka variables are the same 
as in the original model. In contrast, the bka variables 
have been modified to now include the elastic phase shift 
i5cr(e) in Eq. ([5]). This point will be expanded on in 
Sec, mil when we discuss the calculation of the current. 

We turn now to inelastic effects, which arise from 
quasiparticle interactions in the Fermi liquid theory. 
These interactions can be written in a Hamiltonian 
form 3 



tti/'^Tk 



E 



a-<a-' ,{ki} 



^ a<c7'.{k,} i 



X2 



(9) 



where : : denotes normal ordering and v ~ \/ {hvp) is the 
density of state for ID fermions moving along one direc- 
tion. To summarize, the Fermi liquid theory is generated 
by the Hamiltonian H = Ho + i^int, given by Eqs. ^ 
and (HI), with the elastic phase shift ([6]). In fact, Eqs. ([6]) 
and @ correspond to a systematic expansion of the en- 
ergy [3, [sHi, compatible with the SU(N) symmetry and 
the Pauli principle. It includes all first and second order 
terms in the low energy coupling strength cx 1/Tk- 

The great advantage of the Fermi liquid approach is 
that it can also be applied to non-equilibrium situations. 
Note that the Fermi level ep appears twice in the above 
equations: it defines the reference for energies in the elas- 
tic phase shift ^ and also for the normal ordering in 
Eq. ([5]). When the system is put out-of-equilibrium, for 
instance when each lead has its own Fermi level, e p loses 
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(f>i or (f>2 



FIG. 1: Examples of Hartree diagrams for the self-energy built 
from Eq. Q. The full dots (resp. black and grey) indicate 
vertices with four or six external lines, a, (5 and 7 denote 
spins. 



its meaning as a Fermi level and becomes merely an ab- 
solute energy reference. This can be used to relate Q the 
coefficients (ai, a2, ipi, (t>2, X2) as we shall show below. 



C. Kondo floating and perturbation theory 

To make progress in calculating physical observables 
at low energies, we will treat the interacting part iJjnt 
(c.f. Eq. ^) of the Fermi-liquid Hamiltonian perturba- 
tively. Among the various diagrams built from Eq. ([5]), 
it is convenient to separate the trivial Hartree contribu- 
tions to the electron self-energy from the more compli- 
cated diagrams. The former are obtained by keeping an 
incoming and an outgoing line and by closing all other ex- 
ternal lines to form loops as shown Fig. [T] The resulting 
diagrams are then in correspondence with the diagrams 
describing scattering by a local potential. Therefore they 
can be included in the clastic phase shift, 

Tk Tj, \Tk 



2T2 



(10) 



where we have defined SNq^o- — J deSna-{s) and 6Ei^„ — 
J deeSriaie). Sna{e) = na{e) — 9{eF ~ e) is the actual 
quasiparticle distribution (nr,(e-k) = ^\j>kaY) relative to 
the ground state with Fermi energy ep. We see again 
that ep sets the reference in Eq. (fTO)) for both e and 
brifjie). Including Hartree diagrams is essentially tanta- 
mount to a mean-field treatment of the interaction term 
Eq. On a physical level, these Hartree terms can 

be interpreted as a mean-field energy shift of the Kondo 
resonance arising from a finite quasiparticle population 
and their interactions. We shall see that in the case of 
an SU(4) Kondo quantum dot, these Hartree terms play 
a significant role in determining the non-linear conduc- 
tance; this is not the case in the more conventional SU(2) 
Kondo effect. 

While the idea of perturbatively treating i/jnt is 
straightforward enough, a possible weakness of the Fermi 
liquid approach is the number of seemingly undetermined 



parameters in Eqs. ([6]) and ([9]). The standard Fermi liq- 
uid treatment of the Kondo effect allows one to relate 
the coefficients a\ and 01 via the so-called 'floating' of 
the Kondo resonance (to be discussed below); these coef- 
ficients correspond to leading-order Fermi liquid correc- 
tions. However, for transport quantities in the general 
SU(iV) Kondo case, we will see that the remaining coeffi- 
cients, corresponding to higher-order corrections, are also 
important. Luckily, these too can be related to one an- 
other using a novel and powerful extension of the Kondo 
floating recently proposed in Ref. It allows to re- 
late the different phenomenological coefficients of Eqs. © 
and ([9]) ; we describe the basic reasoning involved in what 
follows. 

The Kondo resonance is a many-body phenomenon 
that results from the sharpness of the Fermi sea bound- 
ary [3^ . Physically, conduction electrons build their own 
resonance. The structure of this resonance is therefore 
changing with the conduction electron occupation num- 
bers, as Eq. pU]) shows explicitly. However it can not 
depend on ep, which is a fixed energy reference. This 
idea is implemented by shifting the Fermi level £f by Sep 
while keeping the absolute energy sp + e and the absolute 
occupation numbers ncr(e) fixed in Eq. (jlOp . As a result 
e — > e~Sep and Sna{s) 6ncr{e) + 6{e — Se p) — 9{e) . Im- 
posing invariance of the phase shift leads to the following 
Fermi liquid identities 



ai = (iV-l)0i 



N - 1 



a2 



^2 = {N~2)x2, 



(11a) 
(lib) 



where the first relation (jllap was initially derived for the 
general SU(A^) case by Nozieres and Blandin [l|. 

Note that an alternative way to derive Eqs. ([TT|) is to 
insist that the entire structure of the Kondo resonance 
simply translates in energy when we dope the system 
with quasiparticles in a way that corresponds to a sim- 
ple increase of the Fermi energy Q. Nozieres' original 
derivation of Eq. (|llap in the SU(2) case [s^ also used 
this idea, but restricted attention to an initial state with 
no quasiparticles. Eqs. (|llb[) follow when we apply the 
same reasoning to an initial state having some finite num- 
ber of quasiparticles. Note that for SU(2), or a half-filled 
dot (to N/2), a2 = from Eq. Q so that (j)2 = and 
X2 = 0. The ncxt-to-leading order corrections all vanish 
in agreement wi th p revious works on the ordinary SU(2) 
case IllSliO,!!!!. 

It is worth mentioning that the second generation of 
Fermi liquid terms (02, 02, X2) can also be derived in the 
framework of conformal field theory. In Ref. 0] , a single 
cubic Casimir operator is given, which reproduces the 
three terms corresponding to the coefficients a2 , 02 , and 
X2- The identities (jllb[) are then automatically satisfied. 

The floating of the Kondo resonance (and resulting 
conditions) also has an important consequence for cal- 
culations of observables in the presence of a voltage: the 
results will not depend on where one decided to place the 
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(a) 



©1 




FIG. 2: Diagrammatic construction for the independence of 
observables in £f- Crosses correspond to elastic scattering. 
Two- and three-particle interactions are represented by, re- 
spectively, black and grey full circles. Many diagrams in the 
perturbative expansion in Hint & exhibit a dependence in 
ef- Nevertheless, it is possible to gather and combine those 
diagrams to produce e_F-invariant forms. The combination (a) 
that appears in the irreducible self-energy does not depend on 
ef as a result of Eq. (|lla|) . Combination (b) is a second in- 
variant, thanks to Eqs. (|llbp . contributing to the irreducible 
self-energy, (a) and (b) together imply the phase shift invari- 
ance discussed in the text. Finally, the four-particle vertices 
of (c) can always be combined to cancel the dependence in ef 
thanks to Eq. (jllb|l . Apart from (a), (b) and (c), all other 
diagrammatic parts involve energy differences in which the 
reference ef naturally disappears. The combinations (a), (b) 
and (c) can be understood as emerging from Ward identities 
related to the U(l) gauge symmetry. For example, Eq. (|llap 
has been shown [43, |43[ to derive from a Ward identity with 
a vanishing charge susceptibility. 



a vanishing contribution. In contrast to these simplifica- 
tions, 5Eiu in Eq. (jlOp remains generally different from 
zero due to the energy dependence of the 02 vertex in 
Eq. ®. 

On may wonder whether the physical argument of the 
floating of the Kondo resonance, as presented in Ref. Q 
and repeated in this paper, is sufficient to extend the 
results of this paper to higher orders Fermi liquid cor- 
rections. Applying the floating argument to the next 
(third) order, one obtains an incomplete set of relations 
between the coefficients such that some of them remain 
undetermined. In the language of conformal field theory, 
it means that more than one operator is involved at each 
(higher) order. How to relate the coefficients of those op- 
erators is a rather difficult problem. In the SU(2) case, 
a solution was given by Lesage and Saleur [44| . 

We finally turn to the discussion of the Fermi liquid 
model renormalization. Treated naively, the model leads 
to divergences in physical quantities. It is regularized [Slj 
by introducing an energy cutoff D (different from the 
original band width of the model) larger than typical en- 
ergies of the problem but smaller than Tk ■ Energies in 
Eq.® are therefore restricted to the window [— 
The dependence of observables in D is then removed by 
adding counterterms in the Hamiltonian. It is strictly 
equivalent to the introduction of cutoff D dependence in 
the coupling constants [i^ (a^, 0i, etc). The correspond- 
ing counterterms are discussed in Appendix [XI 



III. CURRENT CALCULATION 



dot Fermi energy ep within the energy window defined by 
the chemical potentials of the leads. On a technical level, 
this is because, by virtue of Eqs. any shift Sep of 

the dot Fermi energy will be completely compensated by 
a corresponding shift in the Hartree contributions aris- 
ing from the quasiparticle interactions. This invariance 
is explained in detail in Fig. [2l Note also that this in- 
variance has physical consequences as well: it implies, for 
example, that the current is not affected by the capacitive 
coupling to the leads (in the Kondo limit). 

Given the above invariance, it is convenient for calcu- 
lations to choose the Fermi level such that 

6NQ,a = j deSn^ie) = 0, (12) 

so that any closed fermionic loop built from an energy- 
independent vertex vanishes. For this choice of position, 
6No_^ vanishes which greatly simplifies the phase shift 
expression (|10p . Moreover, the X2 vertex in Eq. ([9|) does 
not contribute to the current and the noise when the 
perturbative calculation is stopped at second order. The 
reason for that is that the X2 vertex is already second or- 
der and can only appear once. Its six legs are connected 
to at most two current vertices so that at least two of 
these legs must connect to form a closed loop implying 



We now outline the calculation of the current using 
the Fermi liquid theory described in previous sections. 
Again, the complete Hamiltonian is H = Hq + iJjnt 
(c.f. Eqs. (|8|9|) ). corresponding to respectively to elas- 
tic and inelastic scattering; the approach will be to treat 
Hint as a perturbation. Slightly abusing terminology, we 
will include all Hartree contributions arising in pertur- 
bation theory in the free Hamiltonian Hq; Hq will thus 
correspond to the elastic phase shift given in Eq. (|10p . 
Contributions to the current which only involve Hq (thus 
defined) will be referred to as the 'elastic current'. Hint 
is then added perturbatively, without Hartree diagrams, 
in order to compute the corrections due to inelastic scat- 
tering. 

A. The current operator 

The current operator at x is generally given by 

^(^) = ^ E {^U^)d.M^) - d,i'Ux)M^)) (13) 

where m is the electron mass. Various expressions can be 
obtained for the current depending on which basis it is 
expanded. It is convenient [24| in our case to choose the 
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basis of scattering states that includes completely elastic 
(and Hartree terms) scattering, i.e. the phase shift PU)) . 
and that correspond to eigenstates of the single-particle 
scattering matrix. Such states will have waves incident 
from both the left and ri ght leads. This is in contrast to 
another standard choice [2al, which is to use scattering 
states which either have an incident wave from the left 
lead, or from the right. We refer to such states as the 
'left/right' states. 

We first discuss our scattering states in first quantiza- 
tion. Eigenfunctions corresponding to the aka- variables 
do not see the dot or the Kondo effect. Using Eq. 
they read 



'sin6l(e'(*=^+''')=" - 
cos6'(e*('=^+'='^ 



X < 0, 



-z{kp+k)x-^ a; > 0, 



(14) 

where — tt/A measures the asymmetry of the coupling 
to the leads, see Eq. ([2]), and the eigenenergies hvpk are 
measured from the Fermi level ep- The situation is more 
complicated for the b^a variables. The associated eigen- 
functions at small a;, close to the dot, depend on the 
complex ground state wavefunction of the Kondo prob- 
lem. They are not known, and in fact it can not even 
be reduced to a one-particle problem. However, we can 
write the eigenfunctions far from the dot. 



cos 6»(e*('''^+'=)"^ - 5fce-'(''^+'=)^) a; < 0, 
sin6'(e-*('=^+'=)^ - 5fce*('=^+'='^) x > 0, 



(15) 

where the S matrix is related to the phase shift pTH) , Sk = 
g2i5(efc) eigenenergy £k = hvpk. The eigenstates 
and (fTS)) have the same energy. They can be combined to 
give the left and right scattering states with the energy- 
dependent transmission T{e) = sin^(20) sin^(5(£)). In 
the SU(2) case (or generally particlc-holc symmetric 
case), Sq — 7r/2 and the system is closed to unitarity 
for symmetric leads coupling. 

We come back to second quantization and project 
the electron operator ipa{x) over the eigenstates 
and (Uni). Conservation of the current implies that I{x) 
does not depend on x. We choose an arbitrary a; < 
far from the dot, II is the current at x and In at —x. 
If / denotes the conserved current, I = II = Ir- The 
combination sin^ 9 II + cos^ 9 Ir leads to the compact 
expression 



1^ sin 26* \a}^(x)ha(x) — a\{—x)Shcr{—x) 



2 cos 29 



\a\(x)aa{x) - at (-a-)a<,(-.T)]), 



(16) 



with K{x) = Efc&fc<Te"=" and SKix) = Ek^kbkae'''^ . 
Physically, operators taken at x {—x) correspond to in- 
coming (outgoing) states [i^ . The second line in Eq. (fTBl) 
turns out not to contribute to the mean current, the noise 
or any moment of the current. 



Before proceeding with the calculation, it is worth not- 
ing that in the SU(2) case, the proximity to the unitary 
situation allows a simpler treatment [13]. The current is 
written / = /„ — Ibs with /„ = quantum or 

thermal fluctuations are included in the backscattcring 
current Ibs which can be written in terms of aka- and 
bkcr operators [HI, [4^. However, the range of application 
of this approach is restricted to the SU(2) case with a 
completely symmetric leads coupling. In any other sit- 
uations neglecting fluctuations in lu is incorrect [l^, HH] 
and Eq. becomes necessary. 



B. Elastic contribution to the current 

We are now in a position to compute the mean value 
of the current in an out-of-equilibrium situation. A 
dc bias is applied between the two electrodes imposing 
fJ-L ~ Mr = eV. Left and right scattering states, corre- 
sponding to CL^ka and CR^ka operators, are in thermal 
equilibrium with chemical potentials /iL and fij^ . Hence, 
using Eq. we obtain the populations 



{bibk'] 
{alak'[ 



= SkM' [cos^ Ohisk) + sin2 9fR{ek)] , (17a) 

= SkM' [sin^ e/L(efc) + cos^ efn{sk)] , (17b) 
5k k' sin 26*, 



{blak' 



/L/R(e) = /(e - A'l/r), 



-ihiek) - fRisk)], (17c) 
(17d) 



for all spins a. Eq. (|12p . that implies a vanishing 
Hartree diagram, is satisfied with /xl = sin^ 9 eV and 
/iR = — cos^ 9 eV. f{e) = (1 -I- e^'^)~^ is the Fermi distri- 
bution. 

The average current is obtained from Eg. pB|) and re- 
produces the Landauer-Biittiker formula [46| 



Ne 



with the transmission 



&T(e)[/L(e)-/R(e)], (18) 



T{e) = sin^ 26 siT?{5{e)), 



(19) 



and the phase shift 



a2 



Tl\ 3 



[eVf sin^ 29 



5{e) = (5o + — £ - 

(20) 

where we have used the identity plbp . a2 = (iV— 1) 02/4. 
Here, the phase shift (5(e) has an extra (V, T) dependence 
due to mean- field (Hartree) interaction contributions (c/. 
Eq. pO|) ). Within the heuristic resonant-level picture, we 
can interpret this as the voltage inducing a quasiparticle 
population, whose interactions in turn yield a mean-field 
upward energy shift of the Kondo resonance. Note that 
the relevant interactions here are not the leading-order 
Fermi liquid interactions described by tpi^ but rather the 
next-leading-ordcr interaction described by 02- 
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At zero temperature, the current can be expanded to 
second order in eV/Tx. The asymmetry and the zero- 
energy transmission are characterized by 



C = cos 29, To = sin^ (5o 



(21) 




with C = in the symmetric case. The current takes the 
form 



(1 - C^)Ne^Vlh 



= Tq — C sin 2So a\ 



eV 



K 



FIG. 3: Second order diagram describing the interaction cor- 
rection to the current from Hamiltonian Eq. ((9}. The open 
circle represents a current vertex while filled black dots corre- 
spond to interaction vertices, a, /3 are spin degrees of freedom. 
The self-energy term is formed by the three lines connecting 
the two interaction vertices. 



eV 
TV 



cos2(5o(l + iC^)T7: - sin2,5o(l - SC^)^ 
12 6 



(22) 



Kcldysh space given in momentum-energy space by 
1 . / Fo Fo + l' 



£k 



Fo — 1 Fq 



S{e - Sk) 



(24a) 



Inelastic contribution to the current 



The Keldysh framework [49[ is well-suited to estimate 
interaction corrections ([9]) to the current. The mean cur- 
rent takes the form 



/ = {Tj{t)e 



-j.Jcdt'H,,,^(t')x 



(23) 



where the Keldysh contour C runs along the forward time 
direction on the branch 77 = + followed by a backward 
evolution on the branch rj — ~. is the correspond- 
ing time ordering operator. Time evolution of I(t) and 
-ffint(i) is in the interaction representation with the un- 
perturbed Hamitonian Hq Mean values (...) are 
also taken with respect to Hq ([8]) with bias voltage, see 
Eqs. (ini). Note that the time t in Eq. is arbitrary 
for our steady-state situation. Finally, in order to main- 
tain the original order of operators in I{t), we take left 
(creation) operators on the J] = — branch and right (an- 
nihilation) one on the 77 = + branch. 

A perturbative study of Eq. (^5]) is possible by ex- 
pansion in ifint and use of Wick's theorem. This leads 
to usual diagrammatics where one should keep track 
of the Keldysh branch index. The lowest order recov- 
ers the results of Sec. IIII Bl describing elastic scatter- 
ing. The next first order gives only Hartree terms al- 
ready included in Eq. (|22p . Hint gives rise in general 
to three vertices with coefficients 0i, 02 and X2 where 
the last two are already second order in 1/Tk- Thus 
it is consistent to keep only 0i in the second order ex- 
pansion in Hint- A typical Green's function is defined 
by GTh'^^x - x',t - t') = -i{TMx,t,m)bHx',t',r,2)). 
For clarity, spin indices are omitted here and below since 
all nonintcracting Green's functions are spin diagonal. 
Nonintcracting Green's functions are 2x2 matrices in 



Gaaik,e) 



Fo Fo + 1 



(24b) 



Gabik, e) = Gbaik, e) ^ 2i7r(aj^&fc) Q ^ d{e - Sk), 



(24c) 



with the Pauli matrix = 



1 




\yFoiek) = 2{blbk)- 



1, and Foisk) = 2(aj.afe) — 1, as given Eqs. pT]) . We wish 
to compute the second order correction from Eq. (|23|1 . 
It involves the self-energy contribution shown Fig [3] and 
defined by 

S^-^n^i - i2) - C*(^i,ii-i2)x 

kiMM (25) 

xG'^r'ik2,h^h)GZ-''ih,ti-t2) 

The causality identity, for t 7^ 0, 

I]++(t) + S— (t) = S]+-(t) + E-+(i), (26) 

is derived by writing the explicit time dependence in 
Eq. It leads to various cancellations, in particular 

for terms where the lines external to the self-energy ([25]) 
bear no 771/2 dependence. The lines that join the current 
vertex to the self-energy in Figl3] travels from x (or —x) 
to (the dot) and the opposite. Thus, using the Green's 
bmction (|24ap in real space (with a = ±1) 

G'^^ (ax, e) = zTTT/e--/^- (^Fo (e) + | " I ^ 1 ) 

(27) 

and the identity one shows that the terms with 

operators taken at x in Eq. p6p give a vanishing contri- 
bution to the current. This is merely a consequence of 
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causality: interaction, which takes place at a; = 0, can 
only affect outgoing current and not the incoming part. 
We are left with the current correction 



A^(A^- l)esin26' / 



')1 ,V2 



(28) 



de 
2^ 



{^Sg+^"' e)!]"^ '"^ (e)Cy;',^ '- {x,e) + c.c). 



The summation over rji and ?/2 gives two terms: (i) one 
includes the combination — S . It gives a contri- 
bution proportional to D exactly cancelled by a countert- 
erm. Details are given in Appendix [A] (ii) the second 
term involves the combination E'' — S ^ and remains 
finite in the limit D +oo. It reads 



Sin 



N{N - C2)e7r / 



2h 



iS + S*) J g(S-+-S+-)(£)z7ri.A/(e), 



(29) 



with C given by Eq. ^ and A/(£) = fi^ie) - Me). 

We proceed further and restrict ourselves to the zero- 
temperature case. The left and right Fermi step functions 
are introduced by going to frequency space for Eq. (j25p . 
and then by using Eqs. (|24ap and p7)) . The result in- 
volves a sum of terms with products of cos^ 6 and sin^ 9. 
Two distinct integrals, 

Ji= f de f de' T"^ de"fi^{e + e' -e") (30a) 

•'Ml •'Mr •^ 

J2= f de f de' T"^ de" h{e + e' - e") (30b) 

•'MR •'MR •' 

corresponding respectively to one- and two-particles 
transfer, appear with the following combination 

cos^ sin^ iJ. ^2JMJ.- J^i^-^') + ^Ml + C^)^ 

4 

(31) 

With Ji = {eVf/6 and J2 = 4(6^)^/3, we obtain the 
current correction 



(1 - C^)Ne^V/h 



cos25o(iV-l) 



01 ey 

Tk 



Wint, _ ... ^(|)^ eVV ( ^^(^ 

12 4 
(32) 

This result can be given a quite simple physical in- 
terpretation along the line of Ref. [2^. The ^\ term in 
the interaction part of the Hamiltonian ^ can be de- 
composed on the left /right operators basis using Eq. ([2]). 
It then describes processes where 0, 1 or 2 electrons are 
transfered from one scattering state to the other. Us- 
ing Fermi's golden rule and cos^ Q sin^ Q + cos^ B sin^ = 
(1 — C"')/8, the total rate of one-electron transfer is eval- 
uated to be 2ri(l - C"') where 



,eV 6\ f eV 



(33) 



From cos* 9 sin* = (1 — C^)^/16, the total rate for two- 
electron transfer is r2(l - where T2 = For 
one- and two-electron transfers, e cos 25q and 2e cos are 
interpreted as the corresponding charge transfered be- 
tween leads [2^. Writing the current correction as 

5/i„t = (ecos2,5o)2ri(l - C*) + (2ecos5o) y(1 - C^f, 
we recover Eq. ([5^ . 



D. Current for SU(2) and SU(4) 



The results of Sees. IIII B[ IIII CI can be extended to 
finite temperature as explained in Appendix|B] We detail 
results for the total current I = Ic\ + 5I\nt in the {N = 
2, m ~ 1) case and {N = 4, to = 1, 2) cases. 

For SU(2), a single electron is trapped on the dot, ai = 
and 02 = 0. The current takes the form 



1 - 



ai 
TV 



(34) 



where /„i = {2e^V/h){l~C^). In the particle-hole SU(4) 
symmetric case with two electrons, ai ~ 3(/>i and a2 = 0. 
The current reads 



/ = Ir 



Tk 



9 



6 



9 



(35) 

where 1^ = {'ie^V/h){l - C^). 

Turning now to the SU(4) case with one electron on 
the dot, one finds that the inelastic contribution to the 
current vanishes identically (c.f. Eq. as the 'effec- 

tive charges' associated with interaction-induced scatter- 
ing events are proportional to cos 2do and hence identi- 
cally zero [2^. The only contribution is thus from the 
elastic channel (c.f. Eq. ([22]) ). yielding: 



aiCeV 
Tk 



a2 



eV 
TV 



(1-3(7^) 



(36) 



where 1^ = {2e'^V/h){l — C^). There is no temperature 
correction up to this order of the low energy expansion. 
The case with three electrons (to = 3) and SU(4) sym- 
metry is related to the one-electron case by particle-hole 
symmetry. The Kondo resonance is thus changed from 
above to below the Fermi energy. The result for the cur- 
rent is then the same as Eq. (j36p . but with an opposite 
sign for the asymmetry {9 — > 7r/2 — 9, C ^ ^C), i.e. 
the roles of left (L) and right (R) leads are exchanged for 
hole transport. 

The differential conductance G{V) = ^iV) obtained 



from Eq. (j36p gives an asymmetric curve whenever C ^ 0. 
Consider the first the strongly asymmetric case, where 
\C\ becomes sizeable. In this case, the asymmetric lin- 
ear cV/Tk correction in Eq. (j36p dominates even at low 
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bias voltage. For strong asymmetry |C| 1, the con- 
ductance measures the density of states of the Kondo 
resonance [sO] at ±eV. The asymmetric hnear term 
thus follows the side of the Kondo resonance and reveals 
that the resonance peak is located away from the Fermi 
level This behaviour is in fact generic to the SU(N) 
case when the occupation of the dot is away from half- 
filling. In the SU(2) case or generally for a half-filled dot 
(to = N/2), the resonance peak is located at the Fermi 
level which suppresses the asymmetric linear term, see 
Eqs. (p4l) and ([SS]) . 

Turning now to the case of a symmetric dot-lead cou- 
pling (C = 0), we see that as expected, the differential 
conductance G{V) is symmetric in V at all dot fillings; 
hence, it exhibits a quadratic behaviour at low bias. In 
the SU(4) case, the conductance obtained from Eq. pG]) 
is predicted to be maximum at y = 0, in agreement with 
results obtained from slave boson mean field theory [23] . 
Within the Fermi liquid approach, and for one electron on 
the dot, this behaviour is at first glance rather puzzling. 
As we have already indicated, in the SU(4) case, the con- 
ductance is completely due to the elastic transport chan- 
nel. Using the heuristic picture provide by the resonant 
level picture (i.e. elastic scattering due to a Lorentzian 
Kondo resonance sitting above the Fermi energy), one 
would expect that the differential conductance should 
increase with increasing voltage, due to the positive cur- 
vature of the expected (Lorentzian) transmission coeffi- 
cient. This picture is in fact incorrect, as it neglects the 
important Hartree contributions discussed in Sec. IIIIBI 
Heuristically, as the voltage is increased, quasiparticle 
interactions lead to a mean-field upward energy shift of 
the position of the Kondo resonance. Because of the rela- 
tion 02 = (4/3)a2, this energy-shift effect dominates, and 
causes the conductance to decrease; without this mean- 
field energy shift, the conductance would indeed exhibit 
a quadratic increase at small voltages. Note that an in- 
correct upturn in the conductance was reported in pre- 
vious works: Ref. [l^ neglected the higher-order Fermi 
liquid interaction parameter (jjo and the resulting mean- 
field energy shift, while Ref. [2J| treated it incorrectly 
(corrected in 0|). Note also that the results for the 
conductance presented in Ref. [l^ only apply to a sys- 
tem with a strongly asymmetric dot-lead coupling. 



with the current fluctuation Aj(t) = I{t) — (i(t)), see 
Eq. (jl6p for the current operator expression. 

Insight can be gained by first examining the strong 
coupling fixed point at zero temperature, with eV <^ Tk 
so that i5(e) ~ Sq. Quantum expectations in Eq. ([37]) are 
evaluated with the free Hamiltonian ([H]) . The shot noise, 

^0 ^ ^^'^''^' ro(l - C^)[l - To(l - C^)l (38) 

is pure partition noise like a coherent scatterer [i^ . This 
result implies a vanishing noise in the particle-hole sym- 
metric case, like standard SU(2), with symmetric leads 
coupling (To = 1 and C = 0). In this specific case, 
the shot noise is only determined by the vicinity of the 
Kondo strong coupling fixed point, that is by the inelastic 
Hamiltonian ([9]) and the corrections to (5o in the elastic 
phase shift The shot noise is therefore highly non- 
linear with S ^ a.i low bias voltage. Since the corre- 
sponding current is close to unitarity, an effective charge 
e* = (5/3) e has been extracted from the ratio of the noise 
to the backscattering current [2^. e* ^ e should how- 
ever not be confused with a fractional charge. It emerges 
as an average charge during additional and indepen- 
dent Poissonian processes involving one and two charges 
transfer as shown by the calculation of the full couting 
statistics [i^. Nevertheless, this charge e* = (5/3) e 
is universal and characterizes the vicinity of the Kondo 
strong coupling fixed point. It can be seen as an out-of- 
equilibrium equivalent of the Wilson ratio. 

In asymmetric situations (Tq 7^ 1 or C ^ 0), the linear 
part (|38[) of the noise does not vanish and even domi- 
nates at low bias voltage. For instance in the SU(4) case, 
To = 1/2 so that To(l-To) = 1/4. This property is quite 
relevant for experiments and may be used to discriminate 
SU(2) and SU(4) symmetries for which the current gives 
essentially the same answer [l^]. In a way similar to the 
symmetric SU(2) case, we can define an effective charge 
from the ratio of the non-linear parts (~ V^) in the noise 
and the current [H, • This is however less straightfor- 
ward to measure experimentally since it requires a proper 
subtraction of the linear terms. 



IV. CURRENT NOISE 

Fluctuations in the current are almost as important as 
the current itself. In particular, the shot noise (at zero 
temperature) carries information about charge transfer 
in the mesoscopic system. The purpose of this section 
is to detail the calculation of the zero-frequency current 
noise, 

S = 2 j dt{Ai{t)Ai{0)), (37) 



A. Elastic contribution to the noise 

Inserting the current operator (fT6|) in Eq. ([37]) . the 
elastic Hamiltonian ([8|) gives a gaussian measure which 
allows to use Wick's theorem, and thus Eqs. ^T7\ . Like for 
the current, we obtain a Landauer-Biittiker formula [i^ 
for the noise with the same transmission (|19p and phase 
shift ([20]) . At zero temperature, it reads 

S=-— / deT{s)[l^T{s)]. (39) 
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An expansion to second order in cV/Tk yields the elastic 
(non-linear) correction to the noise 



a 3 7 



6S, 



{l-C^)2Ne^\V\/h 
with coefficients, 



5s 



(1)6^ 

Tk 



Ss 



(2) 



eV 



Cai sin 2^0 



[l-2To{l-C% 



(2) 



a 



12 

a2 



i(l + 3C2)(cos4(5o + 2 sin(5o sin 3SqC ) 



(40) 

(41a) 
(41b) 



^(l-3C^)sin25o[l-2ro(l-C^)], 



and the total elastic noise reads Sd = So + SSd- The 
first order correction (|41ap gives an asymmetric part to 
the noise for C ^ 0. In a way similar to the current case, 
particle- hole transformation (6q — > tt — (Sg, 012 —0^2) 
reverts the sign of the asymmetry (|41a[l which indicates 
that the Kondo resonance is centered off the Fermi level. 



Inelastic contribution to the noise 



coax) 




FIG. 4: Diagrams for the noise appearing in the first and 
second order expansions in the inelastic Hamiltonian Eq. ([9|. 
Diagrams (c) and (d) give identical contributions. For dia- 
gram (f), the three- lines bubble can alternatively dress the 
bottom Green's function, a, j3, 7 denote spins with a 7^ /3 
and /? 7^ 7. Open dots represent current vertices while filled 
black dots correspond to interaction vertices. 



We follow the same procedure as for the interaction 
correction to the current established in Sec. IIII CI The 
mean value in Eq. (jST)) is taken within the Keldysh frame- 
work, similar to Eq. The correct ordering of / op- 
erators is maintained by choosing time on the rj = + 
branch and time t on the f] = — branch. The perturba- 
tive study of the noise involves diagrams with two current 
vertices instead of one in Sec lIII Cl The resulting calcu- 
lations are therefore similar to those for the current but 
are much more involved on the technical side. The dia- 
grams relevant for the noise at first and second order in 
1/Tk are shown FigU) Noninteracting Green's functions 
are still given by Eqs. 

Three vertices can be built from the interaction Hamil- 
tonian Hint ® with coefficients (f>i, (f>2 and X2- The X2 
vertex has six legs and appears at most once at order 
1/T|-. Topology therefore imposes that two legs among 
the six must connect to form a closed loop. The cor- 
responding energy integral vanishes thanks to Eq. ((T^). 
Apart from Hartree terms already included in the elas- 
tic part, see Sec lIV Al the expansion to first order in 
Hint gives the single diagram Hja. Involving both <f>i 
and (f>2, the interaction vertex in diagram 2) a is char- 
acterized by the energy dependent coefficient 0(£,e') = 
4>i + (/'2(e + e')/22\'. The corresponding noise correction 
is given by 



dede' 



(42) 



^"'^ (£)(/)(£, £')^'"'(e'), 



where A''^-''^ is the building block defined in Appendix [Cl 
and shown in FiglHl We replace A^'^'^'^ by its expres- 
sion (jC9[) (with sin^((5(e)) instead of Tq) and perform the 
1] summation to get 



"'-'int 



(iV-1) 



dede'Af{e) 



(l-C2)22iVe2//i Tk 
X Fo(e) sm'^{d{e))(j){e, e')Af{e') sm2d{e'). 



(43) 



This general expression can finally be expanded order by 
order in eV/Tx- After energy integration, we obtain at 
zero temperature a first and a second order noise term, 



"'-'int 



{l-C'^)^2Ne^\V\/h 
with coefficients. 



6s 



>,i) 



eV 
TV 



Ss 



(q,2) 



eVV 



54tt^^ = C(/)i(iV - 1) sin2(5o sin^ <5o. 



5s\ 



(a,2) 



:-C2(7V-l)[02Tosin2,5o 
6ai(/)iro(l-4To/3)]/2. 



(44) 

(45a) 
(45b) 



Note that these two terms vanish identically for symmet- 
ric leads coupling (C = 0). It can be checked again that 
the first order correction is odd with respect to particle- 
hole symmetry while the second order is even. 

The expansion to second order in Hmt yields the dia- 
grams Ulb-f with two interaction vertices. To be consis- 
tent with the rest of the perturbative calculation, only 
(/>! is kept in each interaction vertex. The contributions 
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corresponding to the diagrams [4lb-f are all calculated in 
Appendix [Uj Finally, the total noise reads 

S = So + SS,, + 5S^^,+dSl, + 2dS^,,,+SS^^,+dSL, (46) 

where the different terms are respectively given by 
Eqs. dSHl), (gOl), dBl), (ICTOl) . (|CT2l) . and (IC6l) . 



C. Noise for SU(2) and SU(4) 

We have also extended the noise calculation to finite 
temperature along the lines of Appendix [Bl In the asym- 
metric case, the results are too cumbersome to be writ- 
ten here. In the symmetric case, the noise was calculated 
in Ref.0] where it was emphasized that corrections are 
rapidly sizeable at finite temperature. Hence the shot 
noise regime is expected only at very low temperature. 
Keeping a zero temperature, we specialize here to the 
experimentally relevant SU(2) and SU(4) cases with one 
electron on the dot, m = 1. 

In the SU(2) case, the noise correction to Eq. ([38|) reads 



to one lead and only weakly to the other (c.f. Eq. ([21)). 
Transport in this limit corresponds to an incoherent tun- 
neling regime where the hopping from the weakly-coupled 
lead to the dot is the limiting process. It can be checked 
from the Eqs. dS]), (lUTOl) . (ICT2ll . (IClBl) . (IC6l) for 

the noise, and Eqs. ([22]) . ([32]) for the current, that the 
Fano factor F = 1 to leading order in 1 — |C|. This 
is of course expected since the tunneling regime gives 
Poissonian statistics for charge transfer. Note that this 
unity ratio holds order by order for the eV/Tn and 
(cV/Tk)^ correction separately. In addition, we also have 
5o/2e/o = 1 to leading order in 1 — |C|. 

In the opposite limit of a symmetric dot-lead couping 
(i.e. C = 0), coherent effects are important to transport, 
and charge transport is generally not Poissonian [ist . 
Note also that in the symmetric case, the non-linear parts 
of both the current and current noise are cx V'^ . We find 
that the generalized Fano factor is given by 



F 



l + sin2(2(5o) + 



9-13sin^(2i5o) 



N-1 



^sin45o 



jV+4 
N-1 



cos2(5n - 2%sin2(5o 



(51) 



6S 



(l-C2)2iVe3|F|//i 



eV 
TV 



5 4 



(47) 



where the terms cancel each other unexpectedly. 

In the SU(4) case, the noise correction has linear and 
quadratic contributions. 



SS 



(1 



C^)2Ne3\V\/h 

2 



eV 
TV 



Cai 



(1 - 2(7^) 



eV 
TV 



7C4)-^(1-3C2) 
6 



(48) 



V. RESULTS AND CONCLUSION 
A. Main Results 

Following Refs [H, [2^ , we define a generalized Fano 
factor F which describes the relation between the non- 
linear current and current noise: 



2e 6F 



(49) 



It is defined as the ratio between the non-linear parts of 
the noise 6S = S — So (c.f. Eq. ([38])) and of the current 
SI = I — Iq, where 



To = (1 - 



2 NToe^V 



(50) 



is the linear current (for eV <^ Tk)- We focus on the non- 
linear noise and current, as it is these quantities which are 
sensitive to the contribution of Fermi liquid interactions. 

Consider first the strong asymmetric case |C| — > 1 
(i.e. 6* — > or — > ^/2), where the dot is strongly coupled 



This Fano factor includes the effect of interactions; we 
have used the important equahties in Eqs. pT|) . Note 
that this result has no explicit dependence on V/T^ '■ it is 
thus a universal quantity characterizing the Fermi liquid 
properties of the strong-coupling fixed point; also note 
that F is invariant under a particle-hole transformation, 
where m N — m. We stress that the fact F ^ 1 
in general is due both to the existence of two-particle 
scattering at the fixed point, as well as to the partition 
noise associated with single-particle scattering. We give 
in Table [T] values of F for different N and m. 
For N +00, Eq. ^ leads to 



F 



3 cos 4(5o + 4 cos 2^0 
4 + 2cos25o 



(52) 



Note that in the large N limit, two-particle scattering 
processes become insignificant [5l| for the current (since 
(j>i and (j>2 scale as ~ 1/A^) and the result is consistent 
with the non-interacting resonant level. In this limit, the 
Wilson ratio is in fact just one [H, HI]. However, the ef- 
fect of two-particle scattering processes seems to survive 
in the current noise through the diagram of Fig. [41(b). 
Heuristically, this diagram represents an enhancement of 
the coherent partition noise already present in the ab- 
sence of Fermi liquid interactions. The small interaction 
parameter </)^ ~ is compensated by the spin sum- 

mation with ~ equivalent diagrams. The effect is 
therefore linear in A^, at the same level as elastic terms. 

The expression of Eq. ([52ll can be checked in two lim- 
iting cases. For i5o 0, it gives ~ 1. Again it corre- 
sponds to the tunneling regime since a small phase shift 
Sq implies a weak electronic transmission Tq = sin^ So . 
When particle-hole symmetry is recovered, ~ '''/2, we 
find F = —1. In this limit, the conductance is close to 
unitarity and interactions play no role since the diagram 
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m 


N 


2 


3 


4 


5 


6 


7 


8 


9 


1 


-5/3 


-0.672 


-0.300 


-0.156 


0.003 


0.156 


0.287 


0.393 


2 




-0.672 


-3/2 


-1.256 


-1.031 


-0.855 


-0.679 


-0.503 


3 






-0.300 


-1.256 


-7/5 


-1.326 


-1.254 


-1.173 


4 








-0.156 


-1.031 


-1.326 


-4/3 


-1.313 



TABLE I: Fano factor F, Eq. for various iV and m. 



of Fig. [IJb) gives a vanishing contribution for 5q = 7r/2. 
The situation is therefore similar to the ordinary SU(2) 
case [2^, where one has Poissonian weak backscatter- 
ing events. In our case though, backscattering events are 
elastic and imply the transfer of only one electron so that 
F = -1. 

We finally turn to the general asymmetric case, C 7^ 0, 
where we focus on the SU(2) and SU(4) symmetries with 
m — 1. For SU(2), the generalized Fano factor 



F 



h - 

3^3 ' 



(53) 



is obtained from the ratio of the noise p7)) and cur- 
rent corrections at zero temperature. We stress that 
this simple result ([53]) is exact and is not restricted to 
small values of the asymmetry C. E g. (1551) indeed bridges 
the symmetric result F = —5/3 [23] to the tunneling 
regime, = 1 in the strong asymmetry limit C —^ 1. A 
different asymmetry correction was predicted in Ref . [23l | . 
This discrepancy may come from the fact that the current 
expression used in Ref. j23i] is not valid outside the sym- 
metric case C = (see discussion at the end of Sec. lIII A]) . 

The SU(4) case for arbitrary asymmetry is more com- 
plicated since the generalized Fano factor ([^^ bears a 
eV/Tii dependence. This is because the non-linear cur- 
rent and noise have both linear and quadratic corrections 
in cV/Tk (resp. quadratic and cubic terms in V) and no 
simplification occurs when the ratio is computed (univer- 
sality is however recovered in the symmetric case where 
the linear corrections vanish). We therefore prefer to 
compute directly the ratio of the quadratic corrections 
with the result 

^...^gS^-j| '-fn^r +C^ ,54, 



2e 



1-3C2 



where JS'^^) (J/^^^) denotes the noise (current) correc- 
tion to second order in cV/Tk- Again the ratio ([M]) 
connects the symmetric case {N = 4, to = 1 in Table [l|, 

= ~ -0.300 dll to the tunnel or strongly 

asymmetric regime where F^"^^ = 1. Expanding Eq. (j54p 
in C, we obtain F^^) ~ -0.300 (1 - 8.336*2) ^^lich indi- 
cates an important correction due to the asymmetry of 
the coupling to the leads. 



regime using an elaborate Fermi-liquid approach. 

We have particularly focused on the case = 4 rele- 
vant to carbon nanotube quantum dots. One important 
characteristics of the emergent SU(4) symmetry is the 
sign change of the leading current corrections (i.e. linear 
in bV/Tk) as a function of the bias voltage when progres- 
sively tuning the asymmetry between the dot-lead cou- 
plings. More precisely, for a strong asymmetry, we have 
recovered a positive linear correction which traduces the 
fact that the Kondo resonance is peaked away from the 
Fermi level; in this case, the conductance measures the 
density of states of the Kondo resonance at zteV where 
the sign changes with the weakly coupled lead. For sym- 
metric couplings, we have demonstrated that the linear 
correction now becomes exactly zero and that the cur- 
rent becomes maximum at = due to interactions 
via the Hartree contributions. In addition, the noise ex- 
hibits a non-trivial form due to the interplay between 
coherent shot-noise and noise arising from interaction- 
induced scattering events. In the symmetric case, in- 
teractions result in a universal Fano factor F w —0.300 
at zero temperature. For a finite asymmetry between 
dot-lead couplings, the current and the noise have both 
linear and quadratic corrections in cV/Tk- Focusing ex- 
clusively on the quadratic corrections, we have derived a 
formula for the Fano factor which extrapolates between 
the symmetric result and the strongly asymmetric result 
F = 1, perfectly reproducing the Poissonian statistics for 
charge transfer in the tunneling limit. 

In the context of the standard SU(2) Kondo effect, we 
have obtained a generalized Fano factor F = —5/3-1- 
8(7^/3 at zero temperature which is not restricted to 
small values of the asymmetry. Finally, in the limit of 
large A^, it is certainly relevant to observe that the effect 
of interactions tends to subsist in the current noise. 
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Conclusion 



APPENDIX A: COUNTERTERMS AND MODEL 
RENORMALIZATION 



To summarize, we have provided a thorough analysis 
of the non-equilibrium transport in the SU(A^) Kondo 



The improper self-energy can be calculated to second 
order in e/Tk following Rcfs. [13, Uli- The result is that 
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the dependence on the cutoff D can be removed by adding 
the counterterm 



H, 



1 



c.l 



2'ki'Tk 



k.k' 



= -ai7^ hi - 



(Al) 



It corresponds 



to the Hamiltonian Ho+Hint, Eqs. ([8]), 
to a renorniahzation of ai —f ai + Sai. 

We wiU now show that the second contribution that 
arises from Eq. (|28p . and that we have discarded in 
Sec. IIII C[ produces a term linear in D exactly can- 
celled by the counterterm (lAip . Using the identity 
- (t) = sgn(f)(S-+(t) - it takes 

the form 



SI 



(2) 



iV(A^- 1)(1 - C^)- 



2h \ttv'^Tk 



Sx 



J dtsgn{t){E-+ - E+")(t) iTn^Af{-t)+c.c., 



(A2) 



where A/(t) is the time Fourier transform of A/(e) = 
/l(£) — /rXe)- Inserting the Fourier transform of i 1, 



/ -iFois)±l)e~ 

J~D * 



sinh(7rrt) 



(2 cos^ 6le-'^^* + 2 sin^ 6'e-'^«*) 



2e 



(A3) 



t 



in Eq. (^5]) with Eq. (|24ap . it can be checked that in- 
termediate values oi t ^ 1/T,1/V give a vanishing re- 
sult for Eq. ((X2|) (integrand is odd in t). Ea. (|M)) is 
therefore dominated by small t ^ 1/D. In that limit, 
E+-(0 ~ _ e'^*)^/t^ and £-+(<) = (S+"(t))*. 

A/(t) ~ A/(0) + tA/'(0) is expanded to first order in 
t since the zeroth order gives an odd integrand and a 
vanishing integral. After some straigthforward algebraic 
manipulations, we eventually find the result 



SI 



(2) 



(1 - C^)Ne/h 
where we have used that 



= — sin 2do / de 



Sai£ 



K 



[/l(£)-/r(£)], 
(A4) 



Im 



+ 00 



{du/u^){l 



e™)3 = 31n(3/4). 



When higher orders in t are included in the expansion, 
corrections to Eq. (|A4p are of order 0(1/1?) and com- 
pletely vanish in the universal limit D +oo. In par- 
ticular, the 0(1) contribution vanishes by symmetry. Fi- 
nally, the counterterm (jAl[) gives an elastic contribution 
to the current that can be computed along the lines of 
Sec. IIIIBl The result compensates exactly Eq. (|A4[) . 

A second counterterm is generated by vertex correc- 
tions. In the spirit of the self-energy calculation, the sin- 
gular contributions (i.e. depending on the cutoff D) to 






FIG. 5: corrections to the four-particle vertex. 



the four-particle vertex are determined from the standard 
second order diagrams shown in Fig. [5l and proportional 
to 0^. This strong dependence on D is removed by the 
counterterm 



IIc,2 — 



r<o-',{fei} 

D 4 



(A5) 



(>i = (A^-2)0f — -ln2 



To summarize, the perturbative calculation of observ- 
ables to second order in s/Tk {s is a typical energy, ^bB, 
ksT or eV), from the full Hamiltonian IIo + Hint + IIc,i + 
Hcfii leads to finite and well-defined results in the uni- 
versal limit D +00. 



APPENDIX B: FINITE TEMPERATURE 
CURRENT 

We briefly outline how the current is calculated at 
finite temperature. For the elastic part (fT8)) detailed 
in Sec. IIII B[ we merely need the Fourier transform of 
Af{s) = /L(e) - fnie), 



A/(t) 



ttT 



27r sinh(7rTt) 



(Bl) 



The derivatives of A/(t), taken at t = 0, give access 
to the integrals with the corresponding powers of s in 
Eq. UHl). 

The inelastic part of the current is detailed in Sec. IIII Cl 
and given by Eq. (|29p . Eq. ([^^ is evaluated at finite 
temperature by Fourier transform to real time t. The 
time contour is then shifted by irj in the complex plane, 
with rjD ^ 1 (but T, eV <C I/??) such as to suppress the 
dependence on the cutoff D in Green's functions. From 
Ea. (|M|) . the result is (for x- = 0) 



sin^ Oe-'"''*) 



sinh(7rrt) 



with a similar expression for Gi,i,^{—t). The intermediate 
integral result 



dt [g-+{t)g+-{-t)g-+{t) g+- {t)g-+ i~t) 



3 

(B2) 
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is used to derive the current interaction correction 

2 



ox fAT 1^ "^^i 

— — — = cos 2dn [N — 1 ) 



X 



12 



2(7rr)2 



(B3) 



The t-integral in Eq. (|B2p is obtained by first expanding 
the numerator in powers of e^'^^'^''^. Each term gives an 
integral. The standard method to evaluate such integrals 
is to shift the integration contour by —i/T in the complex 
plane wich encloses the pole at a; = 0. 



In order to simplify the forthcoming expressions, let us 
define the following prefactor 



Sp 



_ hN{N - l)sin2 26' 



f— y 

\2vhJ 



(CI) 



We start by considering the diagram [4lf where the self- 
energy bubble is inserted in the top Green's function. 
Going to energy space and integrating over time t in 
Eq. (j37p . the corresponding contribution takes the form 



Si 



Y^ri,ri,J ^^vur,^,)pm,m^,)^ (C2) 



APPENDIX C: DETAILS ON THE 
INTERACTION CORRECTION TO THE NOISE 

We discuss in this Appendix the interaction correc- 
tions to the noise with two interaction vertices, i.e. the 
diagrams Sib- f. Terms oc 02, X2 in the interaction Hamil- 
tonian ([9]) are already second order in 1/Tk. Therefore 
only the term cx is kept for the diagrams |3)b-f since 
the whole calculation goes up to second order in \/Tk. 



where S'''^'''^ (e) is the self-energy part ([25]) that already 
appeared in the calculation of the current. V'''^'^^ is a 
notation for the product of the three Green's functions 
(of the form Q^^^Q^^'^Q'^ ) that enclose the self-energy 
in diagram [Df. Since the current operator (|16p has four 
different terms, this gives a sum of 16 terms for 'P^^-'^^. 
Yet nine of these terms have no 771/2 dependence and 
vanish when summed over ?7i/2. This is a consequence of 
the causality identity Finally ^''i'''^ reads 



+ {-s)g-,^^{~x,e)gz-'[-^.e)gt-{2x,e) + {-sfg^,^^{-x,e)gzHx,e)Qt-i^r^) 

+ {-S*){-S)g-,'^'{~x,e)gll^{x,e)gt~{Q,e) + {-S*)g^;^' {x,e)gZ^ {x,e)g+- {~2x,e) 
+ (-5)g,7^ (-X, (-^. ^)Qta (2:c, e). 



(C3) 



^bb y •^T'^l^bb 

The noise contribution with a bottom self-energy insertion gives a similar expression. Green's functions are replaced 
by their expression ([24]) and the summation over 771/2 is performed together with the causality identity ([26)) . In 
analogy with the current calculation, two sorts of terms are obtained: (i) those including the combination — 
and (ii) those with E"' or E '". Type (i) terms are dominated by energies on the order of the model cutoff D. They 
are exactly cancelled by the counterterm (|A1[) . We therefore only keep type (ii) terms. Combining top and bottom 
self-energy insertion diagrams, (55/^. = 5S(^:^ + ^^(ia.^ ^'^'^ contribution 



5SL ^{N-l)4>l fds 



^([{1 - C')icos4So - cos2<5o)[A/(e)]2 + 2ro(^^o(£)^o(£) - l)]x 

V (C4) 



{l-C^)2Ne^/h 2 v^Tl 

(E+-(e) - E-+(e)) - {F,{e) - l)E+-(e) - {F,{e) + l)E-+(£)) 
At zero temperature, Fo(£)/'o(e) - 1 = -(1 + C^)A/(e) and 

as we have shown in Sec. IIII Cl for the current. The last two terms in Eq. (jC4p involve the combination 4Ji(l + C"' 
^2(1 — C) with Ji/2 given Eqs. ([30]) . Finally, we obtain for the noise correction (jC4p . 



Hl.-^^{^V)\ (C5) 



6SL 



{l~C^)2Ne^\V\/h 





2 


'1 


^4 


/ 5 C2\ 


if 


) 


4 ^ 


T2 + 


V24 8 J 



iC6i 



Next we turn to the diagram|3)c with the particle-hole pair polarization bubble, 

0"-"^^ = E gZ'''Hh,t)g^,r^ik2,-t). (C7) 
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FIG. 6: Building block appearing in diagrams [Ija-e defined 
as A^^''^'^{£), see Eq. (|C9|l . It is formed by one current vertex 
on the branch 77 = + or — with two lines. The incoming line 
connects the branch 772 to ± while the outgoing connects ± 
to 771. 



Another causality identity, similar to Eq. (|26p. also ap- 
plies here. For t ^ 0, 



u++it) + u—{t) = n+-(i) + u-+{t). 



(C8) 



An elementary building block that appears in dia- 
grams IHa-e is shown Fig [6l It consists in one current 



vertex on the branch 77 = -f or — supplemented by one 
incoming and one outgoing lines. In energy space, it 
reads 

A''-"^(e) ^ {i^vf sin20A/(e) [4ToFo(e) + ^771 - 5*772] , 

(C9) 

where the four terms of the current operator (jl6[) are in- 
cluded. Using the definitions ([CT|) . ([UT]) . ([C9|, the noise 
term due to diagram IHc can be written 



5Sf, 



1)1 ,»?2 



(CIO) 



Following a now familiar pattern, there are terms with 
n++ — n and others with 11^'^. The former ones 
depend linearly on the cutoff D and are exactly cancelled 
by counterterms. This will be discussed at the end of this 
Appendix. 



J 



We are left with 



{l-C^)2Ne'^/h 



7r(jV- 1)(1-C^)0^ f dei f de2 
2^ 



27 ^^/(^^)^/(^^) 



(CIl) 



[(n+", +n7+, )cos4Jo + 2(n+"^ , 

L\ ei— £2 £1— £2' ^ ^ El— £2 £1— £2 



n: 



)ro cos25o(Fo(ei)-^^o(£2))] 



At zero temperature, the second term in the brackets gives a vanishing contribution. Developping in terms of Fermi 
step functions, the integrals over energies can be performed leading to 



(l-C2)27Ve3|y|//i ^^'^^ 
corresponding to the combination J2(l — C*^) + 4JiC^. 



(C12) 



r 



The diagram Hid gives exactly the same contribution, 
SSf^^ = SS^j^^. The calculation for the diagram [He is 
quite similar with the introduction of the particle-particle 
bubble, 

11"-"^^) = J2 S',r'iki,t)g'd''Hk2,t), (C13) 
that satisfies the same causality identity (|C8|) . The noise 



term reads 
6Sr„, = 



Sp mm 



dei f de2 



27r j 27r 

xn"i^"=(£i+e2)A"i^"He2), 

leading to 



(C14) 



J 



and 



{l-C^)2Ne^/h 



(n++,, + II-+,J + 2(n++,^ - n^+jTo cos2<5o (Foiei) + F,{e2)) 



-{N- 1)0? (^)\l-C')(l + ^+Tocos2SoC'], 



{l~C^)2Ne^\V\/h 



(C15) 



(C16) 
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corresponding to the combination J2(l + C^) — 4 JiC^ + 8Tq cos 26o J2 — 2 Ji). We finally consider the diagram[4lb 
that also involves the particle-hole bubble H, 

5St,, ^ -Sp{N~l) J2 mV2 f ^A^^^^He,)n^^^^H0)A^^^^^{e2), (C17) 

with the definitions (jCip . (jC7[) and (|C9|) . Keeping only the terms with 11^'^, we obtain 

- (l-C-)sin-(2.o) ^^^:5^^ r / |^A/(s)V(n- +n-). (C18) 



(1 -C2)2iVe2/;i ^ ^ ^ i/2T|. y 

At zero temperature, n^j"^ = IIq ^ = {tti>'^ /2){eV){l — C^) so that the noise correction for diagramlUb finally reads 

SSL _{N-imsm^{25o) ^f^V _ ^2^2^ (C19) 



{l-C^)2Ne^\V\/h 4 \Tk 

Before concluding this long Appendix, we briefly discuss the remaining terms resulting from the 11++ — 11 and 
n++ — combinations. They lead to contributions that are linear in the cutoff D, with corrections scaling as 
0{1/D) and therefore vanishing in the universal limit. The calculation is straightforward and uses the same ingredients 
as in the Appendix [XI i.e. small t dominate time integrals. Therefore we can use IT"' (t) ~ — i/'^(l — e*^*)^/^^, 
n+-(i) ~ -n+-(t) and n-+{t) = (n+-{t))*, n-+(t) = (n+-(t))* in those integrals. The final result reads 

il-C^r2Ne^\V\/h = -SMN-l)Csin2S, sin^.„-, (C20) 

where we recall that S(f)i = {N — 2)(/)i^^ln2. This contribution is exactly cancelled by the counterterm (|A5[) 
included in the diagram of Fig. [Da. 
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